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Achieving efficient terahertz (THz) generation using compact turn-key sources operating at 
room temperature and modest power levels represents one of the critical challeges that must 
be overcome to reaUze truly practical applications based on THz^~^. Up to now, the most 
efficient approaches to THz generation at room temperature - relying mainly on optical rec- 
tification schemes - require intricate phase-matching set-ups and powerful lasers'*'^. Here 
we show how the unique light-confining properties of triply-resonant photonic resonators 
can be tailored to enable dramatic enhancements of the conversion efficiency of THz gene- 
ration via nonlinear frequency down-conversion processes. We predict that this approach 
can be used to reduce up to three orders of magnitude the pump powers required to reach 
quantum-limited conversion efficiency of THz generation in nonlinear optical material sys- 
tems. Furthermore, we propose a realistic design readily accesible experimentally, both for 
fabrication and demonstration of optimal THz conversion efficiency at sub-W power levels. 
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The possibility of enhancing the power radiated by an electromagnetic current source by em- 
bedding it in a resonant cavity featuring both a high quality factor and small modal volume, has 
been recognized since the early days of modem electromagnetic theory^- ^ as an attractive approach 
to extract efficiently the energy emitted by different classes of electromagnetic sources. With the 
emergence of nonlinear optics, this approach was also adopted as a way to enhance a broad vari- 
ety of nonlinear optical phenomena^' ^. Later on, the advent of concepts such as photonic crystals 
(PhC)^°^^^, along with the rapid development of improved micro- and nano-fabrication techniques, 
enabled the demonstration of wavelength-scale counterparts of traditional nonlinear optical de- 
vices, which offered the possibility of realistic on-chip active control of light with Ught itself 
One of the most remarkable illustrations of this concept resides in on-chip enhancement of the con- 
version efficiency of nonlinear frequency mixing processes whose conversion efficiency is intrin- 
sically very low. Of special interest, due to its considerable importance for applications, is the case 
of efficient terahertz (THz) generation. Recently, different approaches to resonant enhancement 
of difference-frequency nonUnear coupUng processes, including those involving a final frequency 
in the THz regime, have been proposed^^"'^. In this paper, we present a scheme that enables en- 
hancement of THz power generation via second-order nonlinear frequency downconversion by up 
to three orders of magnitude compared to conventional nonresonant approaches. By using a com- 
bination of accurate numerical simulations and a rigorous coupled-mode theory, we show how the 
unique properties of photonic microresonators to confine Ught in subwavelength volumes for many 
optical periods can be tailored to dramatically reduce the pump power at which the theoretical max- 
imum of THz generation predicted by the Manley-Rowe quantum limit^ is reached. Furthermore, 
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we show how a generaUzation of the canonical phase-matching conditions enables efficient THz 
generation in nonlinear optical materials in which conventional phase-matching techniques are dif- 
ficult, or even impossible, to implement. Finally, we propose an implementation of highly-efficient 
compact on-chip THz sources operating at room temperature and pumped by sub-//J energy pulses, 
thus opening the way towards practical realization of a number of important applications in a broad 
spectrum of fields, ranging from communications to biomedical imaging^"^'^^. 



To gain physical insight into the approach introduced in this manuscript, we start by revisiting 
the enhancement of the power radiated by an arbitrary current source Jt (r , t) embedded in a single- 
mode lossless dielectric cavity. We assume that the electric field inside the cavity can be written 
as E(r,t) = ar(t) exp{—iujTt) 'E''p{r)/\/U (where E!j^\r) is the electric field profile of the 
resonator mode, U — J dr eo n|.(r)|Ey^(r)p is the electromagnetic energy stored in the resonant 
mode, and arit) is the corresponding slowly-varying electric field ampUtude^^). If we also assume 
that Jt(i", t) = Jr(r, t) exp(— t) and that all the radiated power is collected by means of a 
waveguide coupled evanescently to the cavity, the emitted THz power (Pt) can be written in the 



undepleted regime as 
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where Qt and Qt,s stand for the total quality factor and the external quality factor (i.e., the one 
governing the decay into the waveguide) of the resonator, respectively. Vr = V^/IAt/wt)^, where 
Vt is the effective modal volume of the resonator; whereas ET,max denotes the maximum value of 
|ET(r)|. 
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From Eq. (4), the enhancement of the power radiated by Jt(i", t) inside the cavity is apparent 
through the factor Qt/Vt. Note that a Qt/Vt enhancement factor for sources in resonant cavities 
is famous from Purcell enhancement of spontaneous emission^'^, but also appears in broad range 
of purely classical contexts, from microwave antennas inside cavities to acoustic effects such as 
a blown flute or bowed violin. As we discuss in detail below, it is precisely this enhancement 
factor, together with the unprecedented large values for Qt/Vt that can be realized in photonic 
microresonators^^'^^, that enable increase of the conversion efficiency of optical nonlinear fre- 
quency conversion processes to an extent that can not be achieved by means of any other currently 
known physical mechanism. Note also that Eq. (4) reflects the fact that in order to maximize 
the output power Pt, the intrinsic absorption rate inside the cavity must be matched to the decay 
rate to the waveguide (i.e., for a given value of the Q-factor characterizing the intrinsic absorption 
losses inside the cavity, Qabs, the product Qt/Qt,s x Qt appearing in the right-hand side of Eq. 
(4) is maximized when Qt,s = Qabs)- This Q-matching condition for the efficient extraction of the 
power generated inside the cavity is a general property that appears in a broad range of different 
contexts (thermal, electromagnetic, mechanical, etc) that require efficient emission of energy from 
a high-(5 system^^'^^'^^. 

In the case of a general nonlinear frequency mixing process involving a pump and idler 
frequencies (denoted hereafter by ui and uj2, respectively) in a x*^^^ crystal, the current distribu- 
tion JT(r, t) arises from the temporal variation of the nonlinear polarization vector JT(r, t) — 
(?P^^(r, t)/dt. In particular, in the case of difference-frequency generation (DFG), if we assume 
that the electric fields at oui and U2 are given by Ei(r, t) — ai{t) exp{—iu}it) E^°^(r)/-\/C7i and 
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E2(r, = a2{t) exp{—iu2t) 'E^\r)/^/U2, respectively (where Ui is the electromagnetic energy 
stored in the system at frequency uji), the expression for the power Pt emitted at the final frequency 
takes the following form 

= — ^^T^r^|ai(t)| |a2(t)| |/3eff| (2) 

where (3^^ represents the nonlinear coupling strength between the electromagnetic fields involved 
in the nonlinear difference-frequency mixing. This magnitude can be written as 

/3eff = ' _L ^ " ' ^ (3) 

^Jdvnliv) |Ef)(r)P ^jdvnUv) |Ef (r)^ 
where Xj-j^(r) stands for the spatial distribution of the second-order nonlinear susceptibility tensor 
(here the subindices {i, j, k} stand for the components {x, y, z}, respectively, of the corresponding 
electric field vectors). 

An important figure of merit for the scheme proposed in this manuscript is the enhancement 
factor (77enh) of the output power Pt predicted by Eq. (5) with respect to value of Pt that one 
would obtain using traditional nonresonant approaches to enhance the conversion efficiency of a 
DFG process. In particular, we find relevant the comparison of our scheme with the case in which 
the nonlinear coupling coefficient (/?eff) is increased simply by reducing the nonlinear interaction 
area of the modes involved in the frequency mixing process, using for instance a waveguide for 
the pump, idler and final frequencies (an approach that has been used extensively in the past for 
enhancing the conversion efficiency of DFG in different wavelength regimes*^' ^). For simplicity, in 
this comparison we assume operation in the undepleted regime (i.e., we take ai{t) — ai(0) and 
02 (i) = a2 (0) in Eq. (5)). Additionally, we assume that in our system both the pump and idler fields 



are temporally confined, this confinement being characterized by Q-factors Qi and Q2 for loi and 
UJ2, respectively. As we discuss below for a particular potential implementation of our scheme, this 
temporal confinement not only permits enhancing Pt, but also allows introducing a generalization 
of the canonical phase-matching condition^, which leads to efficient DFG process even in systems 
in which the implementation of standard phase matching techniques (such as birefringence) is 
difficult or even impossible. Keeping these assumptions in mind, after some algebra, one finds that 
it is possible to write an accurate analytical approximation for r^enh as 

^ 2^AiA2 nTnT,en ni,eff n2,e£f ^wg Qt Qjiq q 
^enh ~ ^2,^^ ^2^2 ^2^ Q^^ 1 2 

where Ai, A2, and At stand for the wavelengths corresponding to the pump, idler, and final fre- 
quencies, respectively, ni, ^2, and tit denote the value of the refractive indexes at ui, U2, and ut, 
respectively, of the dielectric medium where the frequency mixing takes place. The parameters 
JT-i.eff; ?^2,eff> and riT.eff represent the effective refractive indexes of the waveguide configuration for 
the corresponding modes at uj-i, oj^, and lot, respectively. A^g is the transversal area of waveguide 
system, while Lwg stands for the corresponding length of the waveguide. Note also that, in order 
to make a meaningful comparison between the cavity system and its waveguide counterpart, when 
deriving Eq. (7), we have assumed that in both cases there exists an optimal overlap (which fulfills 
the phase-matching condition for each case) between Ei(r, t), E2(r, t), and Er(r, t). 

Equation (7) summarizes well the comparison between the different mechanisms that play a 
role to enhance the conversion efficiency of a DFG process in a waveguide configuration (consider- 
ing the optimal case in which phase-matching condition between the waveguide modes is satisfied) 
and the scheme to enhance the conversion efficiency introduced in this manuscript. In particular. 



we have found that, as we numerically demonstrate below for a particular structure, in the case 
of THz generation via a x*^^) DFG process in realistic nonlinear optical material configurations (in 
which considering the absorption losses at the final THz frequency is a key aspect in determining 
the ultimate conversion efficiency), it is possible to reach values 77enh ~ 10^. We emphasize that this 
large value for r]enh is obtained even if the corresponding phase-matching condition is satisfied in 
the waveguide system, which is often challenging to implement due to the vast difference between 
the wavelengths corresponding to the pump and idler modes and the final THz modes. Thus, we 
beUeve that the scheme analyzed in this work could be of paramount importance for increasing the 
conversion efficiency of nonlinear frequency conversion processes whose conversion efficiency is 
intrinsically low due, for instance, to the lack of a phase-matching mechanism in the considered 
frequency range, or due to the small value of x^^^ of the materials of interest. This conclusion is 
the first important result of this manuscript. 

In order to explore the extent to which this concept can be applied to solve the current lack 
of efficient sources and detectors operating at room temperature in the so-called THz frequency 
gap^'^, we illustrate its implementation in a specific structure based on a triply -resonant nonlinear 
configuration. Figure la displays a schematic of the proposed system. The power carried by two 
NIR beams of wavelengths Ai and A2 (playing the role of idler and pump beams, respectively, 
their corresponding power being Pun and P2m) is coupled, by means of an index-guided waveg- 
uide, to two high-order whispering gallery modes (WGM) supported by a dielectric ring resonator. 
These WGM at Ai and A2 are characterized by angular momenta mi and m2, respectively. The 
ring resonator also acts as a dipole-like defect at At, when embedded in an otherwise perfectly 
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periodic THz-wavelength scale photonic crystal (PhC) formed by a square lattice of dielectric rods 
(see the corresponding electric field profile in Fig. lb). Thus, the x^^^ nonlinear frequency down- 
conversion interaction that takes place between the two NIR WGM's circulating inside the ring 
resonator yields a current distribution that radiates inside the PhC cavity at the frequency differ- 
ence ujt = 1^1 — 1^2', the rate at which the radiation is emitted is strongly enhanced by the PhC 
environment in which the ring resonator is embedded. In order to extract efficiently the THz out- 
put power (Pt) from the PhC cavity, we introduce into the system a PhC waveguide created by 
reducing the radius of a row of rods (see Figs, la and lb). In addition, in order to break the degen- 
eracy existing between the x-and y-oriented dipole defect modes, the radius of two of the nearest 
neighbors rods of the ring resonator is reduced with respect to the radius of other rods in the PhC. 
This configuration permits having a large value for factor (Qt/Vt), along with a high-Q resonant 
confinement also for the pump and idler frequencies. 

Figure lb shows the structure that results from optimizing the geometrical parameters of the 
system for efficient generation at 1 THz, assuming a pump beam of wavelength Ai=1550nm, an 
idler beam with A2=1542nm, and that the structure is implemented in GaAs (in which the relevant 
component of the nonUnear susceptibility tensor is du=274pmfV). Note that because the resonant 
cavity in our design forms a filter for the infrared (IR) input power, extremely tight control of 
closely spaced source-frequency lines is not required to produce a specific THz frequency, as 
long as the input spectrum overlaps with the resonant filter bandwidth (although greater frequency 
control and hence greater overlap leads to higher efficiency). Importantly, we also point out that in 
order to maximize the strength of the nonlinear coupling coefficient that governs the energy transfer 
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between the pump, idler and THz fields, the whole structure must be designed so the dependence of 
the THz electric field profile on the azimuthal coordinate inside the ring resonator (see definition 
of 9 in inset of Fig. lb) cancels the modulation introduced in the nonlinear susceptibility tensor 
by the local variation of the pump and idler fields with respect to the axes of the nonlinear crystal 
29,30_ ^Yiis modulation is given by the dependence on 9 of the product E^°\r) E2°'**(r), which in 
the case of two considered WGMs is given by a factor exp [i9{m2 — mi ± 2)]^^(see inset of Fig. 
Ic). For GaAs, and for the above cited values for Ai and A2, we have found that this condition is 
fulfilled by two WGM with mi=572 and m2=575, and a dipole defect mode in the THz-scale PhC. 

The above discussion can be viewed as a generalization of the canonical phase-matching 
condition often found in nonlinear optics^. Specifically, note that in standard phase-matching tech- 
niques, the dispersion relation for the three frequencies involved in the DFG is exclusively deter- 
mined by the intrinsic properties of the particular material in which the nonlinear interaction takes 
places. Consequently, in these traditional schemes, the overall efficiency of a DFG process relies 
entirely on finding a suitable nonlinear material whose dispersion relation permits fulfilling simul- 
taneously, for the frequency range of interest, both the frequency-matching and the phase-matching 
conditions (or alternatively, on finding some physical mechanism, such as quasiphase-matching, 
that permits to phase-match the different fields involved in the nonlinear process). However, in 
the approach introduced here, the dispersion relation corresponding to the final frequency ut is 
different from that corresponding to ui and UJ2 and, importantly, it can be tailored almost at will 
simply by modifying the geometrical parameters that define the THz-scale PhC. This introduces 
a general and versatile route to phase-matching that does not depend exclusively on the intrinsic 
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properties of naturally existing nonlinear optical materials. This fact could be particularly relevant 
in those systems in which the canonical phase-matching condition can not be fulfilled. 

To compute accurately the nonlinear optical dynamics of the structure sketched in Fig. la, 
we have applied a temporal coupled-mode theory (TCMT) formalism similar to that described in 
Refs. 17,25-27. It has been extensively shown that this theoretical framework permits 
characterizing accurately several nonlinear frequency mixing processes, including those in which 
there exists a large difference between the wavelength corresponding to the pump and the final 
frequency^^. Figure 2 summarizes the results obtained in the continuous- wave (cw) regime. In 
these calculations we have assumed that Pun — -P2m (the dependence of the results on the ratio 
P2in/Piin is discussed below) and quality factors Qi = = 3.5 x 10^ and Qt = 10^. These val- 
ues for Q are compatible with both the absorption coefficient of GaAs at 1 THz (Q;=0.5cm~^)^^ and 
the experimental values for the quality factor obtained in similar configurations for the considered 
ring resonator and the photonic crystal cavity^^'^"^. As shown in Fig. 2a, for values of Pii„ > 0.0 IPq 
(see the conversion efficiency (defined here as ratio between the output power at THz and total 
input power at NIR frequencies) starts departing from the conversion efficiency predicted by the 
undepleted approximation, eventually reaching the maximum value predicted by the Manley-Rowe 
relation ^. Specifically, from the steady-state solution of the TCMT equations, one can find that the 
maximum conversion efficiency can be written as rjmax = [1/(1 + ^1/^2)] {^i/^t)^t^i, where 
Tr = Qt/Qt,s and Fi = Qi/Qi,s (see horizontal gray Une in Fig. 2a). This fact is confirmed by 
our numerical simulations (see Fig. 2a). As clearly shown in Fig. 2a, at the critical value of Pun at 
which this maximum conversion efficiency is reached (in our case Pfj„ = 0.19Po> or equivalently 
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^lin =0.32W) the pump power that is coupled to the ring resonator is completely down-converted 
inside the system to power at THz and idler frequencies, giving rise to a sharp minimum in Pur 
and a maximum in P2tr- In the general case in which Pn^ ^ P2in-, we have found that the criti- 
cal powers at which the maximum conversion efficiency takes places (denoted as Pf^^ and P^ir)-: 
satisfy the following expression 

Po ~ 4u;i 

This function is plotted as a solid red line in Fig. 2b. As clearly demonstrated in Figs. 2a and 2b 
(see vertical dashed line in both figures), the maximum conversion efficiency occurs at the intersec- 
tion of P2in = Plin with the expression given in Eq. (6). Note also that that as Pq oc I/Q1Q2QT, 
one can adjust the value of Pfj„ just by varying the product QiQ2Qt- On the other hand, applying 
an analysis similar to the one used in Ref. 17, we have also found that there exists a large region of 
the parameter space {P2m, Pim} in which the considered structure presents multistable response, 
i.e., where there are more than one steady-state solution for a given value of the ratio P2in/Piin 
(see blue area in Fig. 2b). We also point out that the net effect of the absorption losses in the 
conversion efficiency consists simply in downscaling the results obtained in the lossless case by a 
factor TtTi (see dotted line in Fig. 2b). 

In order to completely characterize the THz generation process in the analyzed structure we 
have also studied the temporal evolution of the response of the system to Gaussian pulse excita- 
tions. In these calculations we assume that the temporal width of the pulses corresponds to the 
Ufetime of the THz-scale cavity (tthz ~ 16ns). The value of tthz is much larger than the lifetime 
of the WGM modes at the pump and idler frequencies 0.8ns). Thus, we expect similar conver- 
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sion maximum effiencies as those found in the above cw analysis. Figures 3a-c show the results 
corresponding to three representative values for the peak power of Pun{t) (labeled as A,B and C, 
respectively, in Fig. 2b). As seen in Fig. 3a, when the maximum Piin{t) is equal to the critical 
power P{^^, the pump input pulse is completely consumed after spending approximately 60ns in 
the system (i.e., Pur ~ after 60ns); and, simultaneously, the power of the transmitted idler pulse 
{P2tr{t)) value reaches a peak value that is approximately twice the peak value corresponding to 
the input idler pulse (P2m(^))- For peak values of Piin{t) much lower than the critical power Pf^, 
the almost undepleted behavior can be clearly observed (see Fig. 3b): the peak powers of the pump 
and idler pulses pulses are barely modified as they travel through the system. On the other hand. 
Fig. 3c clearly shows how, for input NIR peak powers well beyond the critical power, the pump 
pulse is completely consumed at i 40ns (see inset of Fig. 3c). However, in contrast to the case 
displayed in Fig. 3a, Fig. 3b shows how after t fti 40ns the up-conversion process that mix ut with 
to yield coi starts being relevant, and, consequently, Pitr{t) begins increasing with t; which in 
turn reduces the overall THz conversion efficiency of the process. Finally, Fig. 3d displays a sum- 
mary of our time-dependent simulations in terms of the ratio between the output THz energy and 
total input NIR energy (defined as Ethz = JZo PAt), and Enir = dt [Pun{t) + P2m(i)], 
respectively). As displayed in Fig. 3a, the maximum conversion efficiency can be reached for 
an input energy Enir =0.04/xJ, which represents a reduction in Enjr of three orders of magni- 
tude with respect to the most efficient schemes for THz generation in nonUnear crystals reported 
up to date ^. Furthermore, note that, in addition to powerful lasers, current efficient schemes for 
THz generation require intricate phase-matching set-ups, whereas in the system introduced in this 
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manuscript maximum theoretically possible efficiency can be achieved in an integrated structure 
having a total area of approximately Imm^. 

In conclusion, we have shown the dramatic enhancement of the conversion efficiency of 
general difference-frequency downconversion processes enabled by triply-resonant photonic res- 
onators. In addition, we have presented a generalization of the canonical phase-matching condi- 
tion, which makes possible efficient DFG processes even in the case of nonlinear optical materials 
where standard phase-matching techniques are difficult or impossible to implement in certain fre- 
quency regimes. By means of detailed numerical simulations, we have illustrated the relevance of 
the proposed scheme by demostrating complete conversion to THz energy of a 0.04//J NIR pump 
pulse in a realistic Imm^-footprint on-chip structure created from GaAs. Alternatively, we have 
demonstrated that in the continuous-wave regime the pump powers required to reach quantum- 
limited conversion efficiency can be reduced to up three orders of magnitude with respect to the 
conventional approaches for THz generation employed up to date. In contrast to previous high- 
efficiency THz generation schemes, the concept introduced in this manuscript opens, for the first 
time, the way to efficient THz generation from sources that are compact, turn-key, and low-cost, 
which we believe could enable a broader use of THz sources. 
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Figure 1 : (a) Schematic of the triply-resonant nonlinear photonic structure analyzed in the text. 
Piin and P2m denote the input powers at the pump and idler frequencies, respectively; whereas Pur 
and P2tr represent the corresponding transmitted powers through the structure. Pthz stands for the 
THz ouput power, (b) Electric field profile (Ez) as obtained from FDTD calculations corresponding 
to the resonant mode appearing at ITHz in the structure shown in Fig. la. The value of the different 
geometrical parameters displayed in this figure are a=102/im, (ii=40.8/im, d2=5.3nm, ds=A.Ofim, 
and w=0.8/im. Inset displays an enlarged view of the electric field profile (pointing along the z- 
direction) corresponding to a whispering gallery with m=572 circulating inside the dielectric ring 
shown in the main figure. The geometrical parameters defining the ring resonator are also shown 
in the inset, Rext and Rint being 40.1;um and 30.5/im, respectively. Yellow areas in both the main 



and inset figures represent GaAs regions, while white areas represent air. 
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Figure 2: (a) Ratio between the total output power emitted by the system at 1 THz and the total 
input power at the NIR pump and idler frequencies. The results for the three frequencies involved 
in the considered nonlinear down-conversion process are displayed (cui, uj2, and ujt correspond to 
the pump, idler, and final THz frequencies, respectively), (b) Solid line renders the dependence 
between P2in and Pun that yields the maximum THz conversion efficiency in the analyzed con- 
figuration. Blue line represents the case = Pun- Yellow and blue areas represent the regions 
of mono-stability and multi- stability, respectively, in the space of parameters {Pun, P2in}- Labels 
A,B, and C correspond to the time-dependent analysis displayed in Figs. 3a, 3b, and 3c, respec- 
tively. 19 
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Figure 3: Analysis of the temporal response of the system shown in Fig. la to gaussian excitation 
pulses. Panels a, b, and c correspond to the peak values for Pun shown by labels A,B, and C, 
respectively in Fig. 2b. In Panels a, b, and c, the scale on the left vertical axis corresponds to the 
NIR power of the pump and idler frequencies, both for the input and transmitted pulses; whereas 
the scale on the right right vertical axis corresponds to the THz output power. Inset of Fig. 3c, 
shows an enlarged view of the temporal dependence of Pitr{t) between t = and t = 50ns. Panel 
d displays the ratio between the total output energy (Ethz) and the total NIR input energy (E^jr) 
as a function of E^jr. Horizontal line in the inset displays the maximum possible conversion 
efficiency given by the Manley-Rowe quantum limit. 
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